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ABSTRACT 

Context. RR Gem is one of the few Blazhko RR Lyrae that has photometric observations available extended enough to study the 
long-term courses of its pulsation and modulation properties in detail. 

Aims. We investigate the pulsation and modulation properties and the relations between them in RR Gem using photometric observa- 
tions from the past 70 years in order to gain further insight into the nature of the Blazhko modulation. 

Methods. We studied the photographic, photoelectric, and CCD light curves obtained at the Konkoly Observatory and other authors' 
published maxima observations. Detailed analysis of the light curves, maximum brightness, and O - C data are carried out. 
Results. RR Gem showed modulation most of the time it was observed. The modulation amplitude showed strong variations from 
the undetectable level (less than 0.04 mag in maximum brightness) to about 0.20 mag. The amplitudes of the amplitude and phase 
modulations showed parallel changes, thus the total 'power' of the modulation have changed during the past 70 years. Parallel changes 
in the pulsation and modulation periods occur with a dP^^d/'iPpuk = 1-6 ± 0.8 x 10' ratio. We also detected 0.05 - 0.1 mag changes 
in the mean maximum brightness and mean pulsation amplitude. 

Key words. Stars: individual: RR Gem - Stars: variables: RR Lyr - Stars: oscillations - Stars: horizontal-branch - Techniques: 
photometric - 



1. Introduction 

There is no satisfactory and commonly accepted explanation 
yet for the phenomenon of the light curve m odulation of RR 
Lyrae stars, the century's old Blazhko effect ( lBlazhkolll907h . 
Less than a dozen Blazhko RR Lyrae have photometric observa- 
tions available that are extended enough to study their long-term 
modulation properties in detail. RR Gem (a = 07^21™33^53, 
6 = -)-30°52'59'.'5, J2000) is one of these stars that has been ob- 
served at the Konkoly Observatory from time to time since 1935. 
The observations were always obtained with the most up-to- 
date photometric technique available at the time: between 1935 
and 1953 photographic, from 1954 to 1983 photoelectric, and 
in 2004 and 2005 CCD observations. (The photographic, photo- 
electric, and the 2005 CCD observations have not been published 
yet.) 

RR Gem was claimed to be Blazhko-modulated with a mod- 
ulation period of 48 days based on the photograph ic observation s 
taken in the 1930s at the Konkoly Observatory (lBalazslll960l) . 
Although the modulation period was dubious, the fact of modu- 
lation itself seemed to be unambiguous. In con trast, in the 1950s 
the modulation seemed to cease (De tr3IT970l) . We continued to 
observe RR Gem in 2004 and 2005 and found it to be modu- 
lated but with very low amp htude and with the shortest perio d 
(7.2 1 6 d) known at that time jJurcsik et al.ll2005at[S6do3l2006h . 

In the first part of this series of papers, the results of the 
2004 CCD observations of RR Gem were reported dJurcsik et alj 



* Tables 4-15 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-FA/| 



l2005al Paperl hereafter). This second part reports our results 
of the reexamination of all the available photometric data of 
RR Gem in order to reveal the long-term behaviour of the mod- 
ulation. 



2. Data 

We utilise photometric observations made by different instru- 
ments of the Konkoly Observatory. Most of these data are pre- 
viously unpublished archive observations. We also use maxima 
timing observations of the past 1 14 years from the literature. 

2.1. Photometric data 

The light curves we use were obtained at the Konkoly 
Observatory between 1935 and 2005. Only a limited number 
of photoelectric data have been published by o ther authors 
(Fitch et al. 1966; lEpsteinl 1 1 9691: Lstepienl Il972t iLiu & .TanesI 
1989; ESA 1997). These data are, however, few in number and 
too sparse for studying the modulation behaviour 

The log of the Konkoly observations is found in Table [T] 
The photographic, photoelectric, and 2005 CCD data, listed 
in Tables 4-12, are only available electronically at the CDS 
( http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+A/). In each table, 
column 1 lists the HJD of the observations. In Table 4, column 2 
gives the photographic Bpg magnitudes, while in Tables 5-12 
differential U, B, V, Rc, h magnitudes of RR Gem are given with 
respect to GSC 02452-01847 and BD -1-3 1° 1547 ^ for the photo- 
electric and CCD data, respectively. 
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Table 1. Log of the Konkoly observations of RR Gem. 
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2.1.1. Photographic data 

The photographic observations were made with the 16 cm astro- 
graph (/ = 2240 mm) of the Konkoly Observatory at Budapest, 
Svabhegy. Between 1935 and 1953, about 1000 exposures were 
taken on about 80 photographic plates. Unfortunately, many of 
the plates were lost or damaged during the past 70 years. We 
succeeded in recovering 598 measurable exposures on 56 plates. 

The photographic plates were digitized on a Umax 
PowerLook 3000 flatbed transparency scanner with a spatial res- 
olution corresponding to about 1.8 "/pixel. Digital aperture pho- 
tometry was applied to the images using standard IRAFQ pack- 
ages. In this way, we determined the photographic blackenings 
of RR Gem and of the photographic comparison stars. 

A series of comparison stars were chosen from the 20 arcmin 
vicinity of RR Gem to cover its brightness variation range. 
Altogether, 15 appropriate comparison stars were found. We 
measured the instrumental B magnitudes of the comparison stars 
during the course of the 2005 CCD observations. The differential 
magnitudes were determined with respect to the B magnitude of 
BD +31°15470 (B = 10.643 mag) we had derived earUer (see 
Paperl). The Bpg magnitudes of RR Gem were determined by 
a 3"* order polynomial blackening curve fit for the comparison 
stars of each exposure. The most deviating points were omitted 
from the fits. 

To check the consistency between our instrumental CCD B 
photometry and the photographic measurements, we also auto- 
matically computed the magnitudes of the comparison stars for 
each of the exposures by the aforementioned blackening curve- 
fitting process using all the other comparison stars except for 
the most deviating ones. The magnitudes determined in this way 
showed good agreement with the accepted instrumental CCD 
B magnitudes for all the 15 comparison stars. The differences 
have a mean of 0.00 mag with a standard deviation of 0.05 mag. 
Neither colour nor brightness dependency of the residuals was 
found according to our CCD measurements. 

A record of the evaluation of only a part of the photographic 
observations has remained available. It contains 331 photometric 
measurements of RR Gem, 70 of them from damaged or lost 
plates. These 70 points were brought onto our magnitude scale 
with a transformation that was derived from a linear regression 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

This star was used as comparison for the 2004 - 2005 CCD ob- 
servations. In the first paper of this series, the wrong BD number was 
erroneously given (BD +31° 1549) for this star. 



of the 261 common data points. The photographic light curve 
was complemented with these 70 measurements. 

2.1.2. Photoelectric data 

The photoelectric observations were made with different unre- 
frigerated photometers between 1954 and 1983. The photome- 
ters in the Newton focus of the 60 cm Heyde telescope (/ = 
3600 mm) at Budapest, Svabhegy, first employed an RCA 1P21 
and later an EMI 9052 B photomultiplier tube. In the Cassegrain 
focus of the 50 cm telescope (/ = 7500 mm) at the Piszkesteto 
Mountain Station, a photometer was used with an EMI 9058 QB 
multiplier. All the instruments were equipped with conventional 
Schott filters (U - UG2 2mm, B - BG12 1mm + GG13 2mm, V - 
GGl 1 2mm), except for the first observations, which were made 
without filter A standard photoelectric reduction procedure was 
applied. No correction was applied to the differential extinction, 
as the comparison star is as near as 8' to the variable. The colour 
extinction was also neglected because the B - V colour of the 
variable and the comparison match at maximum brightness. The 
filtered measurements were transformed to the standard Johnson 
UBV system according to the regularly determined telescope 
constants from standard star measurements. Measurements in 
the U band were only made on 3 nights in the past few years 
of the photoelectric observations. These data are also published 
electronically but not used in the analysis. For the photoelectric 
observations, the star GSC 02452-01847 was used as compari- 
son. We derived the standard B magnitude of 10.79 + 0.01 mag 
of this star from the 2004 - 2005 CCD measurements. 

2.1.3. CCD data 

The 2004 CCD observations and results have been published in 
Paperl. In order to refine the pulsation and modulation periods 
and to measure B magnitudes of stars in a larger area to calibrate 
the photograph ic data , further CCD observations were obtained 
in 2005 ( S6dorl l2006l) with the same instrumentation (the 60 cm 
Heyde telescope with a 1152x770 Wright CCD detector in its 
Newton focus, which gives a field of view of 24' x 17'). These 
data were reduced and transformed to the standard system in the 
same way as for the 2004 observations (see Paperl for details). 

2.2. TheO-C data 

The changes in pulsation period and the phase modulation prop- 
erties are studied using the 0-C values of the maximum bright- 
ness times and the times of the B = 11.3 mag brightness on the 
rising branch, referred to as (O - C)max and (O - C)m=ii.3, re- 
spectively. The fmax and fm=ii.3 timings of the maxima and rising 
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Fig. 1. The (O - C)max diagram of RR Gem over the past 1 15 years according to the average pulsation period (upper line of points). 
The lower line of points represents the {O - C)n,=ii 3 data of the Konkoly light curves, calculated using the same period. These 
points are shifted downwards by 0.2 d for better visibility (see the axis on the right side). The intervals of Konkoly observations are 
indicated. 



branches were determined from the Konkoly observations. The 
O - C values were calculated according to the corresponding 
ephemerides, which ai^e given in the subsequent sections. 



2.2.1. ((9-C)„,x 

Most of the published time s of maxima were collected in 
the GEOj] database dGrafi' '19221; 'Waterfieldl Il927| iLuizetl 
1 9301: llyanovl 1 1 930a|lbinC hudovichev 1930; Dubiagoi [T93of 
Blazhko 1935; Lange 1935, 1969; Kleissen 1939; Batvrev 1950, 
1962; Mania 1956; Guriev 1958 ; Man del 1960; Fitch et al. 
19661: lEpstein 119691; Tsesevic hI Il969l; Kanischeva & Lange 
197U jLangeetal.' '1976; Braune et alJ ,1977;: Liu & Janes 
T989| JeSA' 1997; Vandenbroere 1997, [19981 11999, 2001, 
2003, 2005; Gensler 1998; Agerer & Huebscher 2000, 2002, 
2003; Huebscher 2000, 2003. ,2005; Agerer etal. .2001 



Huebscher et al. 



2006Pwils et al 



2005 
^2006) 



Le Borgn j |2004l; iLe Borgne et alj|2005 
A few additional maxima times of 



Waterfield (1927) and Graf| (1 19221) were also considered. We 
have left out some outlying visual points where no cause of the 
deviation was found. Altogether, 289 reliable (O - C)n,ax data 
points were used. Further 101 (O-C)jnax points were determined 
from the Konkoly light curves by 3"^ - 5* order polynomial fits 
to the data points around the maxima. The times of the max- 
ima were averaged for the B, V, R^, and bands in the case of 
multicolour observations, since the times of maximum bright- 
ness do not differ significantly in these bands. Times of maxima 
in the U band were not taken into account because of the larger 
uncertainty of these data. The fits also yielded the magnitudes 
of maximum brightnesses, which were utilised in studying the 
modulation behaviour. 



2.2.2. (0-C)„ 



=11.3 



The phase modulation can be studied through the O - C data 
of a particular point of the rising branch or of the maximum 
brightness. Because the time of maximum brightness can be de- 



|http : //dbrr . ast . obs-mip . fr/ | 



termined with less accuracy than the time of a certain magnitude 
phase on the rising branch, we decided to investigate the mod- 
ulation of the 11.3 mag brightness phase on the rising branch in 
the B band. This brightness approximately bisects the range from 
the middle of the rising branch to the mean maximum light. 

The (O - C)m=ii.3 variation measures, in fact, a combination 
of the amplitude and phase modulations. It shows oscillations 
due to the changing slope of the rising branch when the light 
curve change can be described as only a scaling in magnitude 
(exact amplitude modulation), as well as due to exact phase mod- 
ulation. Just to estimate the contribution of the exact amplitude 
modulation in the amplitude of the (O - C)n,=ii.3, test data were 
generated using the mean CCD B light curve of RR Gem and as- 
suming exact amplitude modulation. The modulation properties 
of the test data are compared with the observations of RR Gem 
in Sect.O 

The phase modulation can be studied in the (O - C),n=ii.3 
data only if the effect of the long-term period changes are 
eliminated. Therefore we used different ephemerides for the 
different parts of the observations, as the period changes require. 
The ephemerides applied to construct the (O - C)m=ii.3 data 
for the different parts of the observations are given in the 
corresponding sections. 



The maximum brightness times, magnitudes, and 
fm=ii.3 data are given in Tables 13-15 for the photo- 
graphic, photoelectric, and CCD observations, respectively. 
These tables are only available electronically at the CDS 
(http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-fA/). In all 3 
Tables, columns 1-3 give the fm=ii.3 data, their errors, and 
the times of the maximum light, respectively. In Table 13, 
columns 4 and 5 contain the Bpgmax magnitudes and their errors, 
respectively. In Table 14, columns 4 and 6 list the photoelectric 
ABmax and AVmax magnitudes, and columns 5 and 7 give their 
errors, respectively. In Table 15, columns 4 and 5 list the CCD 
ABniax magnitudes and their errors. 
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3. The O -C diagram 

The variation in the pulsation period during the past 115 years 
can be followed in the O - C diagrams shown in Fig. [T] The 
diagrams were constructed using the ephemerides 

f„,^ = 2412077.521 [HJD] + O?397291066 ■ E, 

f^^ii 3 = 2412077.501 [HJD] + O?397291066 ■ E. 

The period coiTesponds to the average pulsation period of 
RR Gem over the past 115 years. The {O - C)m=ii.3 data cal- 
culated with this single ephemeris shows similar changes to the 
(O - C)n,ax data. 

The (O - C)n]ax diagram indicates steady and abrupt pulsa- 
tion period changes. The period of the pulsation varied between 
the extrema of 0.397253 d and 0.3973 177 d. The transition be- 
tween these values occurred abruptly just in the middle of the 
photographic observations, around JD 2 428 800. After the pe- 
riod jump, th e period was decreasing during the next 20000 days 
(see Fig. 1 in iS6doiil2006l) and seem to remain stable in the past 
ten years. Details of the period change rates and observed period 
values are given in the next sections. 

Note that the shape of the (O - C)n,ax diagram, in par- 
ticular the aforementioned sudden period change, strongly 
resembl es that of the Bla zhko RR Lyrae stars, XZ Cyg 
(Baldw in &Samolvkl l2003h and RW Dra jBalazs & Detrd 
ri962l) . 

4. Light curve analysis 

The Konkoly photometric data were analysed in order to detect 
the modulation if existed, to follow any changes in the pulsation 
and modulation characteristics, and to find possible connections 
between them. 

The observations were divided into several parts whenever 
the different observing methods, the data distribution in time, 
and the changes in the pulsation period allowed. When the phase 
coverage and data quality made it possible, Fourier analysis was 
performed to search for modulation components in the spec- 
tra. The modulation was also examined through the maximum 
brightness and the (O - C)ni=i 1.3 data points calculated using the 
actual mean pulsation values. In this way, the maximum bright- 
ness variation caused by amplitude modulation and the phase 
modulation of the rising branch were studied. 

Data analysis was perfor med using the different applications 
of the MUFRAN package (iKoUathl [l990 ) and the linear and 
nonUnear curve fitting abilities of gnuplofl and Mathematicc0 
(IWolframlll996h . 

4.1. Results from photographic observations 

Because of the sudden and remarkable period change that took 
place at the middle of the photographic observations, this data set 
was divided into two parts separated at JD 2 428 800 (denoted as 
pgl and pg2, respectively). In this way, we can follow whether 
the physical processes that gave rise to the sudden and significant 
change in the pulsation period have any effect on the modulation 
properties. 

The pg2 data set has a much less favourable distribution in 
time than that of the pgl data. There are 9 light maxima and 9 
rising branches from two consecutive seasons in pgl, whereas 

http : //www ■ gnuplo t . info/ 

^ Mathematica is a registered trademark of Wolfram Research Inc. 
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Fig. 2. Konkoly photographic light curves for the pg 1 and pg2 
intervals. Solid lines show the mean light curves that take 8 har- 
monics into account. Grey lines show the mean curve of the other 
data set to help compare the pulsation amplitudes. The times of 
the plotted pg2 light curve were transformed in order to elimi- 
nate the period change. 

we have 10 observed maxima and 15 rising branches from the 
15-year long pg2 interval. 

The pulsation period was determined by Fourier fit to the 
light curves taking 8 harmonics into account. The pg2 data set 
does not cover each phase of the pulsation, since the observa- 
tions were focused on the rising branch and maximum light. The 
gap in the middle of the descending branch was bypassed with 
artificial points based on the CCD B mean light curve, shifted 
and scaled appropriately. These points were used only to stabi- 
lize the Fourier fit so they had no effect on any of the conclu- 
sions. The fm=ii.3 data show that a slight period decrease oc- 
cuiTed during the pg2 interval, therefore a linear period change 
was taken into account. The pulsation periods derived for the 
two parts of the photographic data are 

Ppgi =0.397253 + 5- 10"'' d, 

Ppg2 = ^'pg2{) + Ppg2 ■ (t - fpg2 0) d, 

where 

Ppg2o = 0.3973177 + 2 ■ 10"'' d, 
Ppg2 = -(3.7 + 0.8)- 10"'° d/d, andrpg2o = 2428 953.0[HJD]. 

The light curves of the two intervals folded with the appropri- 
ate pulsation periods are plotted in Fig.|2] The fitted mean light 
curves are also shown in the plots. The plotted pg2 light curve 
was trans formed to con stant period by a HJD transformation (see 
details in lJurcsik et al.l l2001). 

From the fm=ii.3 times (listed in Table 13), (O - C)ni=ii.3 val- 
ues were calculated according to the following ephemerides: 

fm=ii.3pgi = 2 427 854.9870 [HJD] + Ppgi -£pgi, 

fm=ii.3pg2 = 2428 953.0176 [HJD]+Ppg2o-£pg2-7.3-10-"-£2g2. 
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for the two intervals, respectively. 

The light curves shown in Fig.|2]indicate that the mean pulsa- 
tion amplitude had increased for the second interval perceptibly. 
No differences in the plate material and in the data evaluation of 
the pgl and pg2 observations account for these discrepancies. 



4.1.1. Modulation properties of tine Pgl data 
(JD2427 833 -2428 248) 

The rising branch of the pgl light curve shows greater scatter 
around the mean curve than that of the pg2 data. Possible expla- 
nations of this phenomenon are incessant period change or phase 
modulation during the pgl interval. Because the {O - C),n=ii.3 
values change sign more than twice during the pgl interval, 
these changes ca nnot be explained with a uni directional period 
change. Recently IChadid & ChapellieJ (l2006l) have reported an 
irregular behaviour of the Blazhko modulation of RR Lyr, and 
they explain this phenomenon by the shock wave passage across 
the atmosphere. According to IChadid (2000), there are impor- 
tant irregularities in the atmosphere of RR Lyr that occur during 
the rising branch of successive pulsation cycles. This effect can 
also be responsible for the widening of the rising branch in the 
pgl data set of RR Gem. The ir regular behaviour on the ris ing 
branch of RR Lyr pointed out bv IChadid & ChapellieJ (l2006l) is, 
however, connected with the modulation. Therefore, as steady 
period change can be excluded, the larger scatter on the rising 
branch of RR Gem indicates the presence of modulation. 

The modulation of the pgl interval is first investigated in 
the rising branch phase data. These data, calculated from the 
fm=ii.3 times (see electronic Table 13) according to the ephemeris 
in Sect. 14.11 were Fourier-analysed. The plausible modulation 
frequency range of - 0.2 c/d were examined. The Fourier 
amplitude spectrum of the pgl (O - C)m=n.3 data is plotted 
in the top panel of Fig. |3] The most significant peak of the 
spectrum is at 0.1 39 c/d frequency. Prewhitening the data with 
this frequency yields the residual spectrum shown in the bot- 
tom panel of Fig. [3] This spectrum only shows low resid- 
ual noise. The best sine curve fit to the (O - C)m=ii.3 data 
gives /modpgi =0.1389 + 0.0004c/d, P^odpgi = 7.20 + 0.02 d, 
and Aphraodpgi = 0.0094 + 0.0004 d. A denotes peak-to-peak am- 
plitude throughout this paper The (O - C)n,=ii.3 data phased 
with this period and the fitted harmonic curve are shown in the 
top left panel of Fig. H] The order harmonic curve fit the 
data with very small scatter, giving convincing evidence of phase 
modulation. 

Due to the large scatter and significant nightly systematic er- 
rors, the maximum brightness magnitudes from the pgl interval 
show the modulation period with uncertainty. These data can- 
not be used to refine the modulation period valid for this in- 
terval, but by accepting the period of the modulation from the 
(O - C)m=ii.3 data, we can estimate the extent of the amplitude 
modulation. The middle left panel of Fig.|4]shows the maximum 
brightness data folded with the 7.20 d modulation period. A sine 
curve fit with this period gives Aampmodpgi = 0.19 + 0.07 mag 
and Bpgi max = 10.83 + 0.02 mag values for the modulation am- 
plitude and the average maximum brightness, respectively. The 
fact that the maximum brightness values can be phased with the 
period derived from the ((9-C)n,=i 1.3 data strengthens the finding 
of modulation during this interval. 




0.00 0.05 0.10 0.15 0.20 

frequency [c/d] 



Fig. 3. Fourier amplitude spectrum of the (O - C)in=ii.3 data 
of the pgl observations. The insert shows the spectral window 
function. The dotted line shows the 4 cr significance level. An 
arrow points to the peak at / = 0.1389 c/d. Bottom panel shows 
the residual spectrum after prewhitening with this frequency. 



4.1 .2. Modulation properties of the Pg2 data 

(JD 2 428 953 -2 434 457) 

The unfavourable data distribution makes the Fourier amplitude 
spectra of the pg2 maximum brightness and (O - C)m=ii.3 data 
uninterpretable because of the severe aliasing effects and the 
possible change in the modulation period. 

Though the period of the modulation cannot be derived from 
the pg2 data, we are able to estimate the amplitudes of the phase 
and amplitude modulations. In the bottom panels of Fig. |4]the 
measured (O - C)n,=ii.3 versus maximum brightness data of the 
individual pulsation cycles are plotted. These plots are indepen- 
dent of the actual value of the modulation period. Any correla- 
tion or connection between these data is the sign that the mea- 
sured scatters of the (O - C),n=ii.3 and maximum brightness val- 
ues are due to modulation rather than observational inaccuracy 
or any irregular behaviour The bottom panel of the pg2 data 
in Fig. |4] shows that this is indeed the case, since the variations 
in (O - C)n,=ii.3 and in maximum brightness are not indepen- 
dent. The possible extent of the amplitudes of the modulation 
can also be read from this plot: Aampmodpg2B = 0.14mag and 
Aphmodpg2B = 0.003 d. If we assume that the pgl and pg2 data 
have similar noise statistics and that the ratio of the modulation 
amplitudes and the scatters of the data are identical for the two 
data sets, then we can estimate the modulation amplitudes of 
the pg2 data simply from the observed scatter of the maximum 
brightness and (O - C)m=ii.3 values (see Table|2]i. The modula- 
tion amplitudes estimated in this way are in very good agreement 
with those read from Fig. |4] confirming that the derived ampli- 
tudes of the modulation are reliable. 

The average maximum brightness was 10.71 ± 0.05 mag ac- 
cording to the pg2 data, which was 0.12 + 0.07 mag brighter 
than observed during the pgl period. Although the scatter around 
minimum light made the estimation of the mean pulsation am- 
plitude somewhat uncertain, it had increased with about 0.15 - 
0.25 mag from the pgl to the pg2 interval. 
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To check whether the observed change of the mean pulsation 
ampUtude and/or the maximum brightness is real or is only an 
instrumental effect, we compared the mean photographic mag- 
nitudes of each of the comparisons stars separately for the pgl 
and pg2 data. The magnitudes show no systematic variations be- 
tween the two intervals, and the magnitude differences are less 
than 0.04 mag for each star. Consequently, we regard the ob- 
served 0.14 mag increase in the average maximum brightness 
as real, and it is most probably a sign of real change in the mean 
pulsation amplitude as well. 

4.2. Results from photoelectric observations 

The photoelectric observations made without any filter are quite 
inhomogeneous as they were obtained during the test phase of 
the new photometer with several different instrument settings. 
Therefore, this data set is only suitable for determining the f^ax 
data but not for investigating the modulation. 

The tm-dx data of the unfiltered observations, the times, and 
magnitudes of all the observed B and V light maxima, as well as 
the rm=ii.3 timing data of the photoelectric light curves, are listed 
in Table 14. 

The pulsation period changed with a constant rate between 
JD 2 435 000 and JD 2 445 000. The change rate has increased 
with about an order of magnitude in the years after JD 2 445 000. 
The time distribution of the observations are uneven as nearly 
the half of the data points were obtained during the first 2 sea- 
sons (pel data, before JD 2 436 629). Later, the observations 
were very sparse: 1-2 nights per annum (pe2 data). The pe2 
data contain the observations made between JD 2 536 956 and 
JD 2 543 931, as the period change rate increased later. 

The period of the pulsation during the pel and pe2 observa- 
tions was 

(f- fpeo)d, 



^pe — ^peO ^pe 



where 



peO 



Ppeo = 0.3973148 + 2 



10"^ d. 



-(7.4 ± 0.2) 



10"'°d/d, 



andfpeo = 2 536 229.0 [HJD]. 



The period and its change rate were determined by a 2"'^ order 
polynomial fit to the Konkoly photoelectric {O - C)ni=ii.3 data. 

The {O - C)n,=ii.3 values were calculated from the pel and 
pe2 fm=ii.3 times (listed in Table 14) according to the following 
ephemeris: 

fm=ii.3pe = 2436 229.0727 [HJD] + Ppeo • £pe - 1 -47 ■ 10-'° • fi^^. 

The pel B data folded with the average pel pulsation period 
and the fitted mean light curve taking 13 harmonics into account 
are plotted in Fig. |5] 

4.2.1. Modulation properties of the Pel data 

(JD 2 436 229 -2 436 628) 

The light curve in Fig. |5] shows no sign of phase modula- 
tion on the rising branch; however, the scatter in magnitude 
around the light maximum is higher than the uncertainty of 
the observations that indicates variation in amplitude, either 
in a regular or an irregular way. The effect of period change 
is negligible during the pel interval. The Fourier amplitude 
spectra of the pel B and V maximum brightness data are 
shown in Fig. |6] Both spectra show two peaks around the 
expected 0.139c/d frequency value. The frequency that fits 
both the B and V data sets the best is /modpei = 0.1374 + 
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Fig. 5. The pel photoelectric B observations phased with the pul- 
sation period and the fitted mean light curve. 



0.0003 c/d (i.e. P,™dpei = 7.28 + 0.02 d). The distance be- 
tween this peak and its neighbour is 0.003 c/d. The -1 cycle/year 
alias frequency results in a fit with significantly higher rms. 
The pel AB and AV maximum brightness data phased with 
^'modpei = 7.28 d are plotted in Fig. |7] The amplitudes of 
the fitted harmonic curves are Aampmodpei b = 0.054 + 0.006 mag 
and Aampmodpei V - 0.052 + 0.007 mag for the B and V data, re- 
spectively. The average maximum brightness in B is 10.77 + 
0.01 mag. According to the fitted mean light curve, the mean 
pulsation amplitude in B band is 1.61 + 0.02 mag. 
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Fig. 6. Fourier amplitude spectra of the pe 1 B and V maximum 
brightness data. Dotted lines show the 4cr significance level. 
Arrows point to the peaks at / = 0.1374 c/d. Panels marked 
with 'pw' show the residual spectra after prewhitening with this 
frequency. Bottom panel shows the spectral window function. 



4.2.2. IVIodulation properties of tine Pe2 data 
(JD 2 436 956 -2 445 673) 



The sparse data distribution in time and the pulsation period 
change during this interval makes these data hard to interpret. 
Therefore, the pe2 data are treated similarly to those of pg2. 
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Fig. 7. Photoelectric B and V maximum brightnesses of the pel 
data folded with the 7.28 d modulation period. 

Upper limits of the amplitudes of the possible modulation are 
estimated from the plots shown in Fig. |4] and from the compari- 
son of the parameters of the pel and pe2 data listed in Table |2l 
The estimations were made assuming identical cr/A ratios for 
the two photographic data sets and also for the two photoelectric 
data sets. 

The scatter of the maximum brightness magnitudes of the 
pe2 data is more than 3 times greater than the average error of the 
individual points. It suggests that modulation at this time might 
also exist, but with even slightly smaller amplitude than in the 
first part of the photoelectric observations. The (O - C),n=ii.3 
data show no sign of phase modulation. The average maximum 
brightness of the pe2 B light curve was 10.76 + 0.02 mag and 
the mean pulsation amplitude was the same as for the pel B data 
within the error range. 

Table 2. Measured scatter of the maximum brightness and 
(O - C)m=ii.3 data and the measured and estimated (italics) am- 
plitudes of their variations. 



data set 


Cr(6max) 


^ ampmod 


a{0 - C\-n=n.i 


'^phmod 




[mag] 


[mag] 


[d] 


[d] 


Pgl 


0.065 


0.19 


0.0034 


0.0094 


Pg2 


0.052 


0.15 


0.0011 


0.0030 


pel B 


0.022 


0.05 


0.0005 


0.0005 


pe2 B 


0.017 


0.04 


0.0004 


0.0004 



4.3. Results from CCD observations 

A detailed study of the modulation properties of RR Gem dur- 
ing the 2004 observations was given in Paperl. Here we only 
summarise the global parameters of the modulation in order to 
compare them with those observed during the photographic and 
photoelectric observations. The fm=i 1.3 and maximum brightness 
data determined for CCD B light curves of the 2004 and 2005 
seasons are listed in Table 15. The (O - C),n=ii.3 values were 
calculated according to the ephemeris: 

fm=ii.3CCD = 2453019.5404 [HJD] + Pccd • E, 



where the pulsation period is Pccd = 0.3972893 + 3 ■ 10"^ d. 

The (O - C)n,=ii.3 and maximum brightness data of 
B light curves are shown in the right panels of Fig. |4] 
These plots were constructed using the modulation period 
valid for the two seasons: PmodcCD = 7.216 + 0.003 d jSodoj 
l2006h . The amplitude of the modulation in (O - C)m=ii.3 and 
maximum brightness are AphmodcCDB = 0.0018 + 0.0002 d and 
AampmodcCDB = 0.114 + 0.004 mag, respectively (see right pan- 
els of Fig. |4]i. The average maximum B brightness of the CCD 
data was 10.76 + 0.01 mag, and the mean pulsation amplitude 
was 1.62 + 0.01 mag. 

5. Discussion of long-term changes 

The long-term changes in the pulsation and modulation of 
RR Gem can be followed in Figs. |4] and |9] The properties of the 
modulation and pulsation at the studied intervals are summarised 
in Table |3] 

Each of the measured parameters (pulsation and modula- 
tion periods and amplitudes, phase difference between maximum 
brightness and (O - C),n=ii.3 data, and mean maximum bright- 
ness) show much larger changes than their uncertainties allow. 

5. 1 . Changes in tiie modulation properties 

In Fig. |4] the modulation in (O - C)m=ii.3 and in maximum 
brightness, as well as the relation between them, are plotted for 
all the studied intervals. The top and middle panels show the 
(O - C)m=ii.3 and maximum brightness values folded with the 
derived modulation periods of the pgl, pel, and CCD data and 
the fitted harmonic curves. Because the modulation period of the 
pg2 and pe2 data could not be determined due to imperfect data 
sampling, these data are plotted versus the HJD of the measure- 
ments. In the bottom panels each point represents a particular 
pulsation cycle for which both fm=ii.3 and maximum brightness 
could be determined reliably. These plots are independent of the 
period of the modulation, and reflect the amplitude relation and 
phase connection of the amplitude and phase modulation com- 
ponents. The plotted ellipses are the combinations of the fitted 
sine curves of the maximum brightness and (O - C)m=ii.3 data 
for the pe 1 , pg 1 , and CCD data. 

It is not easy to make an exact distinction between ampli- 
tude and phase modulation of RR Lyrae stars, because the shape 
of the light curve (i.e. the relative amplitudes and phase differ- 
ences of the Fourier components) always change more or less 
during the modulation cycle. Even if a node on the rising branch 
exists, a slight oscillation in the times of maxima could happen 
(this is o nly the case o f RR Gem at the CCD observations; see 
Fig. 5 in lJurcsik et al.ll2006 ). Especially problematic is how to 
correctly measure phase modulation. As described in Sect. 12. 2. Tl 
the (O - C)m=ii.3 variation measures a combination of the phase 
and amplitude modulations. For comparison, the bottom panel of 
the CCD data in Fig.|4]line shows the (O - C)ni=i u - maximum 
brightness connection predicted from test data of the exact am- 
plitude modulation (see Sect. 12.2.11 1. The actual slope of the ob- 
served (O - C)ni=ii.3 vs. maximum-brightness plots of the larger 
amplitude modulations are flatter than for the test data. During 
the photoelectric observations, when the amplitude modulation 
was the smallest, if there was any, no modulation in ((9-C)m=ii.3 
was detected. 

The bottom panels in Fig. |4]show that significant changes in 
the modulation behaviour of RR Gem occurred during the past 
70 years. The most significant change in the nature of the mod- 
ulation happened at around JD 2 428 800 in conjunction with the 
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sudden change in the pulsation period. Prior to this event, a pro- 
nounced phase modulation was accompanied by the modulation 
in amplitude. After the sudden period increase the character of 
the modulation had changed. The phase of the middle of the ris- 
ing branch show a small variation in excess, if any, to what is 
expected from exact amplitude modulation, while the amplitude 
of the amplitude modulation changes in the range of an order of 
magnitude. The bottom panels of Fig. |4]show that, after the sud- 
den pulsation period change, only the strength of the modulation 
shows significant changes, but its character mostly remains the 
same. 

The parameters plotted in Fig. |9] show that the changes in 
the amplitudes of the (O - C)m=ii.3 and maximum-brightness 
variations are the most tightly connected values. The correlated 
variations of these amplitudes point to the total 'power' of the 
modulation changing with time. Only a few Blazhko stars have 
extended enough photometric data to follow the changes in the 
modulation properties. The best-studied case is RR Lyr itself, 
where the amplitu de of the m odulation shows cyclic changes on 
a 4-year timescale (Szeidl 1976). In other Blazhko stars XZ Cyg, 
AR Her (Almar 1961), AH Cam (Smith etal. 1999), RW Dra 
jBalazs & Detre] ri962). there are definite observations showing 
that the character of their modulation can change considerably, 
although no detailed study comprising thek entire photometric 
data has been published. 

5.2. Changes in the pulsation properties 

The most important and surprising finding of this study of long- 
term changes in RR Gem is the definite change in the mean 
maximum brightness and pulsation amplitude in the middle of 
the photographic observations. We have detected a difference of 
about 0. 12 mag between the mean maximum brightnesses of the 
pgl and pg2 data, while the uncertainties allow only less than 
0.05 mag, according to the comparison stars photometry (see de- 
tails in Sect. 12.1.1b . 

It cannot be excluded that the changes in these parameters 
between the photographic and photoelectric intervals are also 
real, but because of the defects of the calibration of the Bpg mag- 
nitudes, they most probably arise from the differences between 
the photographic and photoelectric B photometries. 

The e xplanation of the Blazhko phenomenon suggested by 
iDziembo wski & Mizerski (2004) involves energy transfer from 
the fundamental radial to a nonradial mode. It predicts that the 
larger the amplitude of the modulation, the smaller the pulsation 
mean amplitude of Blazhko stars. When the modulation was the 
strongest (pgl data), the pulsation amplitude of RR Gem was 
about 0.1 mag smaller t han at any time later, which agree s with 
the model predictions of IDziembowski & Mizerskil ( |2004|) . 

There are very few long-term Blazhko star observations 
available that are homogeneous, extended, and accurate enough 
to investigate similar variations in the modulation and pulsa- 
tion amplitudes. Extended long-term observations in the same 
band have been summarised and/or publis hed for RS Boo (Nagyl 
11998 ', and references therein), XZ Cyg (LaCluvze et al J 12004 
and references therein), R W Dra dSze idl et al. 200^ and ref- 
erences therein), XZ Dr a ( Jurcsik et al.l l2002i and references 
therein), AR Her (| Almai^4 961[, and references thereiri), RR Lyr 
dKolenberg et alj2 006: Smit het alll2003l;ISzeidl et al J 19971 and 
references therein), and RV UMa (lKovacs| ll995', and references 
therein). Nevertheless, only two of them (RR Lyr and AR Her) 
are appropriate for investigation of modulation amplitude and 
mean maximum-brightness changes. 



lAlmaij (1196 Ih showed that in AR Her the brightest maxima 
remained at the same brightness level (see Fig. 12 in lAlmaii 
1961), although the amplitude of the modulation had changed 
considerably. This behaviour indicates that the mean brightness 
was also fainter in AR Her when the modulation amplitude was 
large, similar to what we have found in RR Gem. It supports the 
correlation of the modulation amplitude and the mean maximum 
brightness of Blazhko variables. 

The behaviour of the modulation of RR Lyr is, however, dif- 
ferent. The mean maximum brightness seems to remain unal- 
tered during the 4-year cycle, whil e the modulation amplitude 
changes considerably (see Fig. 6 in ISzeidlll 1 9761 : l^estonll 1 965l 
Fig. 2). Nonetheless, the cyclic behaviour of the variation in the 
modulation amplitude of this star also differs from the irregular 
changes of RR Gem and AR Her, therefore it can have a different 
physical origin. 

The possible differences between the pg and pe data does 
not allow us to make any conclusions about further changes in 
the mean maximum brightness and pulsation amplitude values 
between these observations. During the pe and CCD observa- 
tions, though the amplitude of the modulation has changed these 
properties of the pulsation remained the same, indicating that the 
above-mentioned correlation is not strictly valid. 



5.3. Connections between the pulsation and modulation 
properties 

The period changes of the pulsation and modulation indicate a 
positive correlation between the periods as shown in Fig. [8] al- 
though they can be compared only on the bases of three epochs 
of data. The data define dfmod/dPpuis = 1.6 + 0.8 x 10^ for 
RR Gem, which completes the list of Bl azhko stars with mea- 
sured period change of both their periods dLaCluyze et al.ll2bo4 
Table 8) with a new item. 

It was shown in lJurcsik et al.l(l2005bl) that the largest possible 
amplitude of the modulation depends on the pulsation frequency 
of Blazhko stars. However, at a given pulsation frequency, the 
modulation has very different amplitudes, from very low values 
to the extent of the possible maximum value. That those Blazhko 
variables with long-term photometric observations show signifi- 
cant changes in their modulation amplitudes proves that the mea- 
sured amplitude of the modulation is just a temporal parameter. 
It may explain why the pulsation frequency gives only a limit for 
the possibly highest value of the modulation amplitude and does 
not determine its actual value. 

Most probably the 'instability' of the modulation is a gen- 
eral property of Blazhko stars, and any plausible explanation of 
the phenomenon should give a reason for the irregularity of the 
observed light curve modulation (see also IChadid & ChapeUien 
l2006b. 

What the modulation properties (amplitude, phase relation 
between the amplitude, phase modulation, etc.) tell us about the 
physics of the phenomenon, if they show considerable changes 
on larger timescales, is unclear. During the time of the photo- 
electric observations, RR Gem showed hardly any modulation. 
This fact warns us that the modulation can only be temporar- 
ily detected in some cases. In contrast to previous expectations, 
the amplitude of the modulation can also be small (lJurcsik et al.l 
2006), below the detectability limit of most of the extended sur- 
veys. Therefore, the most plausible conclusion is that the modu- 
lation is a common, intrinsic property of RR Lyrae stars. 
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Table 3. Modulation and pulsation properties derived from the observations. 



data set 


Ppuis (error) 


^mod (error) 


j^puis (error) 


firaax (error)" 


^phraod(error) 


j4ampmod(error) 


A phase(error) * 




[d] 


[d] 


[mag] 


[mag] 


[d] 


[mag] 


[rad] 


Pgl 


0.397253(5) 


7.20(2) 


1.55(2) 


10.83(2) 


0.0094(4) 


0.19(7) 


-1.8(3) 


Pg2 


0.397317(2)" 




1.77(3) 


10.71(5) 


0.003 


0.15 


0.6(7) 


pe unf. 


0.3973156(35) 














pel B 


0.3973148(2) 


7.28(2) 


1.61(2) 


10.77(1) 


0.0005(4) 


0.05(1) 




pe2 B 


0.3973160(2)" 




1.61(2) 


10.76(2) 


0.0004 


0.04 




CCDB 


0.3972893(3) 


7.216(3) 


1.62(1) 


10.76(1) 


0.0018(2) 


0.114(4) 


0.0(1) 



" The errors of the mean maximum magnitudes were determined in two different ways. When the modulation period was 
found, it corresponds to the zero point uncertainty of the fitted sine curve. For the pg2 and pe2 data, the given errors are 
the scatters of the maximum points, which overestimate the true errors, if this scatter originates partially from modulation. 
* A phase denotes the difference between the phases of the maximum magnitude and (O - C)n,=ii.3 variation. 
" As the pulsation period changed during the pg2 and pe2 observations, the given values correspond to the middle of these 
intervals. Details of the period change are given in Sects. 1431 and l4.2l 
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Fig. 8. Modulation period versus pulsation period. 



6. Conclusions 

The rise in atmospheric shock waves during the pulsation of 
RR Lyra e stars is well-e stablished both observationally and theo- 
reticallv. lPrestonl (1 1 965b proposed that the depth of the formation 
of the main shock wave during rising light shifts periodically in 
the course of the Blazhko cycle. Chadid & Chapellier (2006) re- 
port irregular changes in the atmosphere of RR Lyr during the 
Blazhko cycle. They find that the residual scatter of the radial 
velocity curve is largest when the nonlinear effects are the most 
intense, which is during the shock wave passage across the at- 
mosphere. 

The uncertainty of our photographic data and the un- 
favourable data distribution of both the photographic and pho- 
toelectric data sets do not allow us to distinguish between ob- 
servational errors and irregularities. However, the detected sys- 
tematic variations give clear evidence of the light curve modu- 
lation. Analysis of the 70-year long photometric observations of 
RR Gem has revealed some new and important properties of the 
Blazhko modulation. 

RR Gem is the first Blazhko RR Lyrae star where the am- 
plitude of the pulsation showed definite changes during the 
observation times. An unambiguous amplitude increase hap- 
pened together with a sudden pulsation period increase around 
JD 2 428 800. Before this event, a clear and strong modulation 
was observable with an amplitude of nearly 0.2 mag in maxi- 
mum brightness and 14min in phase on the rising branch. After 
the sudden pulsation period change, the strength of the modu- 
lation decreased, along with an increase in the mean pulsation 
amplitude. The photoelectric observations of 1958 - 1959 show 
clear but weak modulation only in amplitude, with a slightly 
longer Blazhko period. The 2004-2005 CCD observations show 
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Fig. 9. Variation of the pulsation and modulation properties 
of RR Gem. The plotted quantities from top to bottom are 
{O - C)n,ax, pulsation period, pulsation amplitude, average max- 
imum brightness, modulation period, modulation amplitude in 
(O - C)m=ii.3, and modulation amplitude in maximum bright- 



modulation with an amplitude of about 0. 1 mag in B band and 
with very weak phase modulation. 

Cyclic changes in the modulation and pulsation properties 
are known for two Blazhko v ariables, RR Lyr (Szeidl 1976) and 
XZ Dra ( Jur csik et alj |2002|) . Supposing that the changes ob- 
served in RR Gem are also a part of a cyclic behaviour on a 
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century-long timescale, both the pulsation and modulation peri- 
ods are expected to decrease, while the amplitude of the mod- 
ulation is expected to increase in the future. If a sudden jump 
in the pulsation and modulation properties of RR Gem, similar 
to that observed in 1937, would happen in the future that would 
give us a unique opportunity to study the connections of the pul- 
sation and modulation properties of Blazhko variables in detail. 
The extremely short modulation period (7.2 d) makes it easier to 
check the modulation amplitude of RR Gem annually. 
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Fig. 4. The modulation in {O - C)m=ii.3 (top), in maximum brightness (middle) and the relation between them (bottom) for the 
studied light curves. The {O - Qm^n 3 and maximum brightness points are folded with the determined modulation period for the 
pgl, pel, and CCD data, and are plotted versus HJD for the pg2 and pe2 data. The fitted harmonic curves are also shown. In the 
second bottom panels, each point represents a particular pulsation cycle for which both (O - C)m=ii.3 and maximum brightness 
could be determined. These plots are independent of the modulation period and reflect the ampUtude relation and phase connection 
of the amplitude and phase modulation components. The plotted elUpses are the combinations of the fitted harmonic curves. In the 
bottom right panel we show the relation valid for test data of exact ampUtude modulation. 



